EOS at FAIR energies and the role of resonances 
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Abstract. Two microscopic models, UrQMD and QGSM, are used to extract 
the effective equation of state (EOS) of locally equilibrated nuclear matter 
produced in heavy-ion collisions at energies from 11.6 AGeV to 160 AGeV. 
Analysis is performed for the fixed central cubic cell of volume V = 125 fm^ 
and for the expanding cell that followed the growth of the central area with 
uniformly distributed energy. For all reactions the state of local equilibrium is 
nearly approached in both models after a certain relaxation period. The EOS has 
a simple linear dependence P = c^e with 0.12 < < 0.145. Heavy resonances 
are shown to be responsible for deviations of the Cg(T) and c^(/xb) from linear 
behavior. In the T-/ib and T-/ig planes the EOS has also almost linear dependence 
and demonstrates kinks related not to the deconfinement phase transition but to 
inelastic freeze-out in the system. 



1. Introduction 

One of the principle questions of the Compressed Baryon Matter (CBM) experiment 
at GSI FAIR is the equation of state of hot and dense matter produced in heavy-ion 
colhsions at energies about 20 - 40 AGeV [1]. Because the perturbative quantum 
chromodynamics (pQCD) is not apphcable to soft processes with smaU momentum 
transfer, one has to rely on microscopic models that correctly describe many features 
of the collisions at various energies. Two of such models, ultra-relativistic quantum 
molecular dynamics (UrQMD) [2] and quark-gluon string model (QGSM) [3], are used 
to extract the effective EOS of the excited matter in heavy- ion collisions at bombarding 
energies ranging from AGS to SPS. The UrQMD and, to a lesser extent, QGSM 
were already employed for studying the equilibration processes, see [4,5]. Recently 
we modified the analysis by extending it to a non-fixed cell, which should follow 
the expanding area of uniformly distributed energy density [6]. By using both the 
UrQMD and QGSM for studies of the relaxation process in a broad energy range one 
can expect that the model-dependent effects, caused by application of a particular 
event generator, will be significantly reduced. - The models use different mechanisms 
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of string excitation and fragmentation. UrQMD relies on the longitudinal excitation, 
whereas the color exchange scheme is employed in QGSM. The fragmentation functions 
that determine the energy, momentum, and the type of the hadrons produced during 
the string decay are also different. Finally, both models do not use the same tables 
of hadrons, chosen as discrete degrees of freedom. Whereas the UrQMD contains 
55 baryon and 32 meson states together with their antistates, the QGSM takes into 
account octet and decuplet baryons, and nonets of vector and pseudoscalar mesons, as 
well as their antiparticles. Heavy resonances are not included in the current version of 
the QGSM, and this circumstance can be used to elaborate on their role in the EOS. 

Central gold-gold collisions with zero impact parameter 6 = fm were simulated 
at bombarding energies i?iab — 11.6,20,30,40,80 and 160 AGeV, respectively. The 
total energy, the net baryon charge and the net strangeness extracted for a certain 
volume of the reaction, were inserted into a system of nonlinear equations [4] to obtain 
temperature T, baryon chemical potential hb and strangeness chemical potential ns 
of an ideal hadron gas in equilibrium. If the yields and transverse momentum spectra 
of particles obtained in a snapshot of microscopic simulations at time t were close to 
the results of statistical model (SM), the matter in the cell is considered to be in the 
vicinity of equilibrium. Then its equation of state can be derived and studied. Because 
the cell is an open system with instantly changing energy and particle density, the 
verification of the equilibrium conditions is repeated after a time-step of A< = 1 fm/c. 



2. RelELxation to equilibrium and EOS in the cell. 
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Figure 1. Hadron yields in the 
central V = 125 fm'^ cell of central 
Au-I-Au collisions at 40 AGeV in 
microscopic models (histograms) 
and statistical model (symbols). 



Figure 2. Energy spectra 
of hadrons in the central V = 
125 fm^ cell. 



In the standard approach the test-volume was a fixed central cubic cell of 
V — 125 fm"^. The yields of some hadron species are displayed in Fig.[l]for central gold- 



EOS at FAIR energies and the role of resonances 



3 



gold collisions at i?iab — 40 AGeV. The agreement between the results of microscopic 
and statistical model calculations is good after t > 9fm/c. Here the standard criterion 
[yield{mic) — yield{S M)] / error {S M) < 1 is applied. According to model analysis, 
after t « lOfm/c almost all many-body processes going via the formation of strings or 
many-particle decaying resonances are ceased, and one deals mainly with elastic and 
quasi-elastic reactions. The energy spectra dN/AirpEdE calculated microscopically 
are shown in Fig. [2l The Boltzmann fit to particle distributions is presented in 
Fig. [2] as well. Both in UrQMD and in QGSM the energy spectra agree well with 
the exponential form of the Boltzmann distributions. Because the hadronic matter 
in the central cell nearly reaches the state of thermal and chemical equilibrium, 
the macroscopic thermodynamic parameters of the system, such as temperature and 
chemical potentials, become meaningful. 

Isentropic expansion of relativistic fluid is one of the main postulates of Landau 
hydrodynamic theory [7] of multiparticle production. As can be seen in Fig. [3] the 
entropy per baryon ratio is nearly conserved in the equilibrium phase of the expansion 
within the 5% accuracy limit. The entropy densities s obtained for the cell in both 
models are very close to each other, but, because of the difference in net-baryon sector, 
the ratio s/pB in UrQMD is about 15-20% larger than that in QGSM. 
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Figure 3. Entropy per baryon Figure 4. Equation of state: 

in the central cell as a function of microscopic pressure P vs. the 

time t. energy density e. 

Any hydrodynamic model relies on the equation of state, which links the pressure 
of the system to its energy density. Otherwise, the system of hydrodynamic equations 
is incomplete. The corresponding plot with microscopic pressures Pmic(£) is presented 
in Fig. m For both models the shapes of the distributions are very close to linear for 
all energies in question. Thus the EOS has a rather simple form 

P{e) = cle , (1) 
where the sonic velocity in the medium Cg = {dP/de)^/"^ is fully determined by the 
slopes of the distributions P[e). To account for possible deviations from a straight line 
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behavior the slopes of the functions P versus e were averaged over the whole period of 
the equilibrated phase. For the UrQMD calculations the velocity of sound increases 
from 0.13 at Ex^\, = 11.6 AGeV to 0.146 at Ex^\, = 158AGeV. It saturates at cl « 0.15 
for RHIC energies [4] . That corresponds to change of the nuclear compressibility from 
140MeV(AGS) to 200MeV(SPS and RHIC). In QGSM calculations the averaged 
sound velocity is about 0.015 units smaller. Note that due to the averaging over 
time, respectively energy density, these values are lower the maximal values for 
that are reached in the corresponding reactions. Both models indicate that at the 
energy around i?iab = 40 AGeV the slope of the c^('\/s) distribution is changing, and 
the velocity of sound becomes less sensitive to rising bombarding energy. 
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Figure 5. The sound velocity 
in the central cell of volume V = 
125 fm^ as a function of baryon 
chemical potential /ib- 



Figure 6. Temperature de- 
pendence of the sound velocity. 
Dashed line corresponds to calcu- 
lations within Hagcdorn model of 
ideal hadron gas. 



Figure [5] shows the dependence of the on the baryon chemical potential fiB- For 
three bombarding energies, Slab = 20 AGeV, 30 AGeV, and 40 AGeV, the functions 
c^(/iB) are close to each other. In QGSM calculations depends linearly on /ie and 
the slope c^/ is unique for all reactions. In UrQMD the picture is more complex. For 
the late stages of system evolution the slopes of all distributions are also similar, but 
for energies of E'lab > 40 AGeV one sees the rise of the sound velocity at the beginning 
of the equilibration, plateau, and the falloff. This can be taken as indication of the role 
of heavy resonances, because their fraction is presented in the particle spectrum at 
the early period and disappeared completely at the end. These resonances are rare at 
Eiah < 20 AGeV, and distributions c^(/.tB) obtained in both models are quite similar. 

The obtained EOS is soft, because for the ultrarelativistic gas of light particles 
the sonic speed is Cg = 1/ a/S. But the presence of resonances in particle spectrum 
generates the decrease [8] of the Cg. Employing the empirical dependence p{m) cx 
m" , 2 < a < 3 [9], where p{m) dm denotes the number of resonances with masses 
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from m to m + dm, one arrives to the equation of state in the form [8] 

e = (a' + 4)P, (2) 

i.e., \ < c^s < \- This result is reproduced in microscopic models. Note that PHENIX 
collaboration reported the value Cg ~ 0.35 ±0.05 [11], i.e., « 0.12 ±0.3, for Au+Au 
collisions at top RHIC energy y/s = 200AGeV. This value is close to our results and 
also implies rather soft effective EOS. 

Temperature dependence of the sonic speed (T) is depicted in Fig. [S] together 
with the EOS calculated in [10] within the Hagedorn model with /j, = 0. For 
-Blab = SOAGeV and 160AGeV the UrQMD data exhibit a falloff in cl{T) at 
T > 120 MeV in accord with the Hagedorn model. This decrease is assigned to heavy 
resonances, because neither the UrQMD calculations at lower energies nor the QGSM 
calculations without the heavy resonances reveal the negative slope in the equation of 
state Cs(T). Below T = 100 MeV both microscopic models indicate rapid drop of the 
sound velocity that occurs much earlier compared to that of the Hagedorn model. 

In the modified analysis the central cell was further subdivided into the smaller 
ones embedded one into another. If the e of the inner cell is not the same (within 
the 5% limit of accuracy) as the energy density of the outer one, the SM analysis of 
the thermodynamic conditions is performed for the inner cell, otherwise the outer cell 
becomes a new test volume. This permits one to follow the expansion of the area with 
uniformly distributed energy. EOS in the T — /xb plane is shown in Fig. [71 Symbols 
and dashed lines show the evolution of these quantities in a cell of instantly increasing 
volume (T^nit = 0.125 fm"^), whereas dotted (upper plot) and full (both plots) lines 
are related to calculations with the fixed volume V = 125 fm^. The transition to 
equilibrium proceeds quite smoothly if the analysis is performed for the fixed cell. In 
contrast, in the area with uniformly distributed energy the transition is characterized 
by a kink distinctly seen in each of the phase diagrams in both microscopic models. 
The effect, which takes place along the lines of the constant entropy per baryon, is 
caused by the significant reduction of the number of processes going via the formation 
and fragmentation of strings, i.e., chemical freeze-out. The observed phenomenon can 
easily mimic the signature of the QCD phase transition in the r-/iB plane. Evolution 
of strangeness chemical potential with T in the fixed volume and non-fixed volume 
is displayed in Fig. [51 As in Fig. [71 all systems develop kinks in the T{fis) distributions 
precisely at the moment of transition from nonequilibrium to equilibrium phase. Both 
baryon density and strangeness density are decreasing in the test volume, however, the 
baryon chemical potential increases with time, whereas the strangeness one drops. The 
evolution of the /ig and /ie with T proceeds quasilinearly, thus reducing the deviations, 
caused by nonzero chemical potentials, of the functions s{T) and s{T) from the ideal 
gas behavior at = 0. 

3. Conclusions 

In summary, both microscopic models favor the formation of the equilibrated matter 
for a period of about lOfm/c for all reactions in question. During this period the 
matter in the central cell expands with constant entropy per baryon. The equation 
of state can be approximated by a simple linear dependence P = a{y/s)e, where the 
square of the speed of sound = a{y/s) varies from 0.13 (AGS) to 0.15 (SPS) in the 
UrQMD calculations and from 0.11 (AGS) to 0.14 (SPS) in the QGSM ones. 
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Figure 7. Temperature T vs. 
baryon chemical potential /^b- 
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Figure 8. Temperature T vs. 
strangeness chemical potential iig . 



Heavy resonances are responsible for negative slope in Cs(T) at T > 100 Me V in 
accord with the predictions of Hagedorn model of hadron resonance gas. At lower 
temperatures both microscopic models indicate a rapid drop of the sonic speed in 
stark contrast with the Hagedorn model calculations with zero chemical potential. 

Study of the expanding area of isotropically distributed energy reveals that the 
relaxation to equilibrium in this dynamic region proceeds at the same rate as in the 
case of the fixed-size cell. However, here both microscopic models unambiguously 
show the presence of a kink in the T-^jl-q phase diagrams. The higher the collision 
energy, the earlier the kink formation. Its origin is linked to the freeze-out of inelastic 
reactions in the considered area. 
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